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We study the sensitivity of weak lensing by large scale structures to the evolution of dark energy. 
We explore a 2-parameters model of dark energy evolution, inspired by tracking quintessence models. 
To this end, we compute the likelihood of a few representative models with varying and non varying 
equation of states. Based on an earlier work, we show that the evolution of dark energy has a much 
stronger impact on the non-linear structure growth than on the angular diameter distance, which 
makes large scale structure measurements a very efficient probe of this evolution. For the different 
models, we investigate the dark energy parameters degeneracies with the mass power spectrum shape 
r, normalisation as, and with the matter mean density fljvf- This result is a strong motivation for 
performing large scale structure simulations beyond the simple constant dark energy models, in order 
to calibrate the non-linear regime accurately. Prospective for the Canada France Hawaii Telescope 
Legacy Survey (CFHTLS) and Super-Novae Acceleration Probe (SNAP) are given. These results 
complement nicely the cosmic microwave background and Super-Novae constraints. Weak lensing 
is shown to be more sensitive to a variation of the equation of state, whereas CMB and SNIa give 
information on its constant part. 

PACS numbers: 98.08. Cq,98.62.Sb,98.65.Dx 



Dark energy is a generic way to describe the acceler- 
ation of the universe. Within this framework, the accel- 
eration of the expansion measured by the type la Super- 
Novae surveys Q, Q is explained by the contribution of 
a new component. In a FRLW metric, this component 
is described by its equation of state. In this paradigm, 
the cosmological constant is one possible model, among 
others, of dark energy. 

The question of the properties of this dark energy 
remains open. Different observations have been pro- 
posed to evaluate them. Measurement of the distances 
to Super-Novae 0, 01 or of the size of structures on 
the CMB 0, 0] provide informations on how the dark 
energy modifies the relation between cosmological dis- 
tances and redshifts. The evolution of large scale struc- 
tures also probe the properties of dark energy. This 
kind of tests can be built using cluster abundances 
S H EL Ly-a fore st 111 , strong [H EH and weak 
lensingjll \wi Hfl Irl ITafT 

In this article, we will investigate the weak lensing con- 
straints in the case of a varying equation of state. This 
case has been first investigated qualitatively by Benabed 
and Bernardeau Q (hereafter BB01). We will expand 
their results and propose the first quantitative analysis 
of the efficiency of the shear two points function to probe 
varying equation of state dark energy models. The con- 
straints on the dark energy using the linear regime alone 
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are not particularly strong, even using the tomography 
technique [ifj- Here, we take into account the non-linear 
regime of gravitational collapse, which is known to con- 
vey much of the dark energy sensitivity E3 , an d we study 
the degeneracy of the dark energy parameters with the 
matter density Qm, the mass power spectrum shape F 
and normalisation erg. 

As said above, large scale structures are sensitive to the 
expansion rate of the universe, which makes the structure 
growth sensitive to the dark energy content of the Uni- 
verse E3- This sensitivity is due to the fact that when 
the dark energy starts to dominate the energy budget 
of the Universe, the efficiency of the gravitational col- 
lapse is reduced. Hence, the density fluctuations growth 
changes when the dark energy differs from a simple cos- 
mological constant. For models within the current SNIa 
constraints, it corresponds to a slower growth in the lin- 
ear regime. The stronger the variation in the evolution 
of dark energy, the earlier this effect occurs and the more 
suppressed will be the structure growth. Keeping the am- 
plitude of density fluctuations fixed today, this translates 
into an earlier entrance of the fluctuations in the non- 
linear re gim e, leading to more concentrated dark mat- 
ter halosjlfll EfT ] , and therefore stronger lenses 0] . As 
shown in BB01, the shear two points function is sensitive 
to these two effects. It provides an unbiased measure of 
the projected density power spectrum in both the linear 
and non-linear regimes, which is a direct test of the evo- 
lution of large scale structures, and through it, a measure 
of the properties of dark energy. In particular, the scale 
at which the transition between the two regimes occurs 
is tested. This feature is the key to constraint the prop- 
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erties of dark energy. 

In the following, we first review the computation of the 
shear two points function with a dark energy component. 
Then, we propose a simple two parameters model that 
encompass the major features of a generic class of dark 
energy, namely the tracking quintessence models. Next, 
we address the question of the efficiency of future lensing 
surveys to determine the dark energy properties. 



quickly damped by the evolution. Moreover, it is ex- 
pected that they are negligible on scales smaller than the 
horizon at recombination jg, EH S| ■ In the following, we 
only consider such scales and can thus safely ignore these 
fluctuations. 

With these assumptions, the impact of dark energy on 
the evolution of the large structures is completely de- 
scribed by the equation of state [14] . 



THEORY 



B. Shear measurements on distant galaxies 



We establish in this section the theoretical basis of the 
result presented in section |HJ We show here how to com- 
pute the non-linear cosmic shear power spectrum with a 
non-trivial dark energy, and how it is evaluated from the 
data. We discuss the sensitivity of this quantity to the 
evolution of an hypothetical dark energy component. We 
then propose a dark energy equation of state parameter- 
ization suitable for a class of dark energy models. 



A. Cosmology 

Lets assume that the dark energy component is mini- 
mally coupled to the universe, which means that it inter- 
acts with the rest of the Universe via gravity only. The 
expansion of the universe is completely described by the 
Friedmann equations 



8tt 



3M Planck 

47T 



3M, 



Planck 



E^x (i) 

E(^X + 3P X )- (2) 



and the knowledge of an equation of state for each com- 
ponent 



w 



(3) 



The radiation, matter and curvature equation of state are 
fixed. The only unknown quantity here, is the equation of 
state parameter of the dark energy, wq . It usually varies 
between +1 and —1. The case wq = — 1 corresponds to 
a cosmological constant. 

It has been proposed recently that wq can also take 
values lower than —1 [2jJ. This is only possible if the 
dark energy has a negative kinetic energy. Such unusual 
behavior have only been found in very specific models [22j, 
and it is not yet clear if it has any physical meaning. 
Therefore, we keep the conservative prior wq > — 1. 

We also assume that the dark energy does not fluc- 
tuate and thus is not coupled to the fluctuations of the 
matter density. This is of course not true in general: it 
is expected that dark energy fluctuates from one Hubble 
volume to the next, which is expected to leave an imprint 
on the metric fluctuations when the universe expands. 
However, it has been showed that these fluctuations are 



The deviation of light by the gravitational potential 
wells distorts the image of the distant galaxies. This 
shear effect can be used to probe the projected mass dis- 
tribution along the line of-sight (see and references 
therein) , from a measurement of the shape of the lensed 
galaxies. The lensing effect produced by the large scale 
structures is weak, but has already been measured [26]. 

The gravitational lensing effect depends on the second 
order derivatives of the gravitational potential projected 
along the line-of-sight. The convergence k and the shear 
7 describe the distortion of the image of the distant im- 
ages (located at some redshift z s ), by the inhomogeneous 
matter distribution along the line-of-sight. At linear or- 
der, convergence and shear field are related 
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At the same order, the convergence in the direction 8, 
which describes the isotropic change of the image at po- 
sition Xs, is given by: 
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xs , i>(x.-x)v(x)s(v(x)o,x) 



where Xs( z s) is the source radial distance located at red- 
shift z s , and a = 1/(1 + z) is the scale factor. The radial 
distance at redshift z is given by 



X{z) 



dz-. 
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The angular diameter distance V is defined by 



sm(VK X )/VK, 

X.- 

smh{VK X )/VK, 



K > 1 
K = 1 
K < 1 



(7) 



where K is the curvature. As a consequence, the weak 
lensing effect is a direct and unbiased measurement of 
the projected density contrast. 

We focus on the convergence power spectrum P K {t). 
It can be shown [23, US] that it is directly related to the 
3-dimensional mass power spectrum P^d via: 
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where the function describe the lensing efficiency, 

»(x.-x)»(x) 
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When the lensed galaxies are distributed in redshift, the 
observed signal is given by EqEjl integrated among the 
source redshift distribution p s (x)- I n that case, the g(x) 
function becomes: 



g(x) = v( x ) dx'p s (x') 
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(10) 



The source redshift distribution p s (z) is normalised, and 
usually parametrised as 



p s (z) = V 



1 
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The free parameters a, (3 and z s are adjusted to accom- 
modate different survey properties. 

The ellipticity of the galaxies is an unbiased measure of 
the shear, from which we derive the statistical properties 
of the con verg ence field (see a review of the observational 
results in [2fl|). In practice, the statistic of interest is the 
aperture mass variance as function of scale, (Affp(#)) 
[2oT |. also called the M a p statistic. It links the variance 
of the convergence field (which is the field of physical in- 
terest, because proportional to the projected mass den- 
sity) and the shear (which is the observable quantity). 
Its main feature is to provide a natural separation be- 
tween the cosmological signal (which is curl-free) and the 
systematics (contributing to the curl mode). It has al- 
ready been measured on several galaxy surveys The 
M a p statistic at a scale 9 C is defined as the convergence 
smoothed with a compensated filter U(8). Using Eq. J1J, 
it is also given by a properly smoothed shear component 
If- 
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with / e ° d9 9U(9)=0, and where 



Q( 9 ) = m M d'u{9')-u{9). 



(14) 



The tangential shear 7< at a location 6 
(Ocos<p, Osinip) is defined by 
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The choice of U{9) is arbitrary provided it has a zero 
mean. In the following, we will use 



U{9) 




(16) 



For this particular choice, the variance of the convergence 
is expressed in term of the shear power spectrum as 
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The variance of the aperture mass is therefore a broad- 
band estimate of the convergence power spectrum given 
in Eq.®, which can directly be estimated from the 
galaxy shapes. This is the quantity which will be used to 
study the evolution of dark energy in section [H] In or- 
der to make predictions on cosmic shear observable, one 
only needs to compute the convergence power spectrum 
in dark energy. As shown in Eq. (JHJ), it only depends 
on the cosmology through the relation between angular 
distances and redshifts, and through the power spectrum 
of the matter density fluctuations. The source distribu- 
tion p(z) can be determined from the data and does not 
depend on the cosmology. 

A complete discussion on the computations of the weak 
lensing power spectrum with a cosmology with a non- 
trivial dark energy has been done in BB01. We only 
reproduce here the conclusions of this work. 



C. Cosmological distances 

The relation between the cosmological distances and 
redshift is given by eq. ©. The dark energy component 
only leaves an imprint on x( z ) by modifying the acceler- 
ation of the universe obtained by solving Eq. I|1I2J) . 

The modification of the relation distance / redshift 
affects mildly the convergence power spectrum. It can 
be summarized into two simple effects : a normalization 
change and a scale shift (similar to the modification of 
the position of peaks in CMB). 

The lensing efficiency function (Eq. 191 Kill acts as a se- 
lection window. It is maximum for lenses located roughly 
at mid-distance between the observer and the source 
galaxies. This selection effect can be approximated by 
replacing the g{x) 2 term in in Eq. (|HJi by a Dirac func- 
tion 



9 2 (x) ^ 9 2 (x m id) S (x ~ X m id)- 



(18) 
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In this approximation, the normalization change adver- 
tised earlier is driven by the change in the position of the 
maximum of the selection function. 

The scale shift is also easily understood with this ap- 
proximation. This shift comes from the Pzd{ d\ x ) ) term 
in Eq. (HJ. The modification of the maximum of the 
efficiency window selects a different depth for projection. 

In the following, we will show that the matter power 
spectrum can be split into two evolution regimes. At 
large scale, the linear regime is well described by a power 
law. The effect of dark energy can be completely re- 
absorbed into a change in normalization. At small scales 
however, in the non-linear regime, the power law approx- 
imation is no longer a good description of the power spec- 
trum. Like the modification of the peaks in the CMB 
power spectrum, the scale under which is seen the tran- 
sition from linear to non-linear regime will be slightly 
shifted by the above effect. 

D. Power spectrum of matter 

From a practical point of view, there is no need to 
compute the whole power spectrum of the density fluc- 
tuation. Only a narrow range of scale, from a few arc- 
second (galaxy scale) to a few hundreds of arc-minutes 
across the sky is enough. Large scales (> 5 degrees) are 
difficult to access observationally anyway; the weak lens- 
ing correlation amplitude is small, where the residual sys- 
tematics might be a problem, and the surveys capable of 
such measurement are not yet planned. At scales smaller 
than a few arc-seconds, the number of lensed galaxies 
drops, and the noise blows up. 

At redshift one, a few degrees corresponds to a few 
hundreds of Mpc, which is far below the horizon size at 
recombination. As stated in section llAl we can safely as- 
sume that for the scales of interest, a classic CDM power 
spectrum is a good approximation of the power spec- 
trum of the matter density fluctuation. Moreover, for 
these scales, the power spectrum behaves essentially as a 
power law. The remarks made above on the cosmologi- 
cal distances holds, and the effect of the modification of 
the equation of state of dark energy on the cosmological 
distances translate, for this range of scales, in a simple 
normalization shift. 

We have yet to investigate the evolution of the power 
spectrum from recombination until now. The growth of 
structures is modified by the presence of a dynamic dark 
energy component. At the linear order, it is given by the 
well known equation[30] 

D+(t) + 2HD+(t) - ^H 2 n (t) D+(t) = 0. (19) 

In this equation, the matter acts as a source term that 
increases the depth of the potential well and tends to 
increase the density contrast. On the opposite, the ex- 
pansion of the universe acts, via the second term, as a 
friction effect and reduces the efficiency of gravitation to 



increase the density contrast. This term carries all the 
effect of dark energy on the growth of structures. 

At high redshift, the dark energy is completely domi- 
nated by the radiation or the matter density. During the 
radiation era, the growing mode of equation Eq. I|19H can 
be obtained analytically 

D+ocai (20) 

During matter domination, the solution is also well 
known 

D+ oc a. (21) 

These remarks allow us to integrate Eq. Ijl9|l m our 
model; we do not know any other initial condition that 
would allow us to perform the integration and be sure to 
keep the growing solution. 

It is worth noting here that there is no other way to 
compute the growth of structure0,|23|. The well known 
integral solution of Eq. ltl9ll is valid when the universe 
only contains matter, radiation curvature and a cosmo- 
logical constant. It is also easy to convince oneself that 
there is no solution to Eq. ltl9|l that can be integrated 
from today toward the past. Equation l|T9jl generically 
admits a growing and a non-growing solution. Only the 
first one is of cosmological interest, and there is no way to 
build this solution when integrating from the final point 
of the evolution. 

During the expansion, the growth follows the radia- 
tion, and then the matter solution. When dark energy 
gets closer to the energy density of matter, the friction 
term grows compared to the source one. The efficiency 
of gravitational collapse to build up the density pertur- 
bation decreases and the growth of structures is damped. 
For the set of models where the dark energy happens to 
dominate earlier, this reduction of the growth rate is ex- 
perienced at a higher redshift. The exact starting point 
of this damping depends on the evolution properties of 
the dark energy model. Models with wq > — 1 experi- 
ence this effect earlier than for wq — — 1 . For a constant 
equation of state, the energy density of the dark energy 
goes as 

Q Q oca- 3(1+t0£i) . (22) 

As said above, the wq = — 1 model is the cosmologi- 
cal constant case. When wq > —1, Qq grows as the 
scale factor a goes to unity. In this case, the dark energy 
contribution to the expansion is significant at a higher 
redshift than when wq — — 1. This is even more impor- 
tant for varying wq, as shown figure The modelling 
of the dark energy used in this figure will be described 
later. 

Figure El shows the result of a numerical integration of 
eq. l(T9jl for different models. The damping of the growth 
appears earlier in the varying equation of state models, 
compared to the constant equation of state. 

For a fixed redshift, the modification of the growth of 
structure at linear order will be degenerated with the nor- 
malization of the power spectrum. Of course, following 
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Figure 1: The energy density of dark energy normalized to the 
critical density as a function of redshift. Thick plain line is 
the classic A model. The thick long dashed line and thin dot 
dashed line are resp. a wq = —0.8 and wq = —0.6 models. 
The short dashed and dotted line are resp. wo = —0.8, wi = 
0.2 and w — —0.8, w 1 — 0.3 models (co mple te description 
of the parameterization can be found sec. II Ell . The sooner 
the dark energy gets close to one, the sooner it will affect the 
expansion and the growth of structure. As expected, models 
with a equation of state different from wq = — 1 contribute 
significantly to the acceleration sooner. Models with a varying 
equation of state contribute yet sooner. A constant wq = 
—0.6 model interpolates between the two wi / ones. 
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Figure 2: linear growth for different models. Models are the 
same than fig. [jThe growth are normalized to the recombi- 
nation era. The modification of the equation of state induce 
a precocious acceleration that decrease the efficiency of gravi- 
tational collapse at higher redshift than in the wq = — 1 case. 
A variation in the equation of state (wi 7^ 0) amplify this 
effect. As expected from fig. [I]a constant equation of state 
model can partly mimic a varying equation of state: if one 
knows the CMB normalization and today normalization of 
the fluctuation of structures, one cannot distinguish between 
a wq = —0.6 and at»( = —0.8, wi = 0.3 models. 



the linear power spectrum along the line of sight allows 
to break this degeneracy. Unfortunately, the measured 
shear power spectrum is only a projection of the mass 
power spectrum along the line of sight. We showed in 
previous section that this projection acts like a selection 
effect on a single redshift plane. Therefore, the integrated 
growth effect will be, in the linear part of the lensing 
power spectrum, indistinguishable from a normalization 
shift. 

At small scales, the evolution of the density contrast 
changes dramatically. Virialized objects are formed and 
evolve in a different regime than the simple one described 
by the linear approximation. BB01 showed that this 
regime can potentially solve part of the degeneracy de- 
scribed above. The weak lensing measurement on galax- 
ies allows access to the bigger scales of this non-linear 
regime. 

The perturbation approach cannot describe this regime 
as the density contrast is very big at the scale of virialized 
objects. A complete computation of this regime cannot 
be done analytically. One has to rely on hypotheses on 
the properties of this regime to be able to describe it. 
Several "classic" description of this reg ime has been pro- 
posed (among others see [3lL 1321 |3jj|). We will follow 
here the choices made in BB01. We assume that the 
Stable Clustering Ansatz provide a valid description of 
the smaller scales of the non-linear regime. It states that 
virialized objects are stable, that is to say that their phys- 
ical size does not vary with the expansion of the universe. 
Hence, at the scale of these objects, the density contrast 
has to grow to match exactly the expansion. Instead of a 
growth of order a or smaller, the density contrast evolves 
as a 3 / 2 . 

One should note that the scales described by the 
Ansatz are below the shear measurement scales. The 
transition between the linear and non-linear regimes is 
described by a mapping between the two regimes pul ]. 
This mapping is calibrated using a fit to n-body simu- 
lations, as described in Peacock & Dodds [32]. At large 
scale, the mapping keeps unchanged the linear regime, 
and at small scales, it must go as (a 2 / g 2 P) 3 / 2 . Although 
this type of mapping has been widely tested for many 
different cosmologies, it has never been tested for dark 
energy models. However, given the robustness of the 
mapping for very different cosmologies, we assumed it 
remains valid for the class of models studied here. 

This is a very strong assumption. It can partly be 
justified by the fact that it is unlikely that a smoothed 
dark energy component with no coupling can affect the 
small scale behavior of the matter. Its influence should 
only appear as a change in the expansion and thus, as we 
have shown above, as a modification of the linear growth 
of structures. The strength of this argument argument is 
enhanced by a recent result from another description of 
the non-linear regime usually called the halo model. This 
approach describes the virialized object as dark matter 
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halos of known 1 profile and abundance depending on the 
cosmological parameters (23- The results and concepts 
behind this approach have been successfully tested in the 
context of dark energy (lfll l2Cj. In particular the differ- 
ences observed between halos in A cosmologies and in 
cosmology with non trivial dark energy can be explained 
by an earlier entrance into the non-linear regime. The ob- 
served discrepancies are, as expected, all explained by the 
modification of the linear growth of structures (14L Il7l| . 

As we just emphasized above, the effect of a non trivial 
dark energy on the non-linear regime are all encoded in 
the linear growth. As the structures grow slowly, they 
reach the non-linear regime earlier (for the same final 
normalization). Thus, changing their growth at a higher 
redshift they undergo a longer non-linear evolution. This 
translates, in the context of Stable Clustering Ansatz into 
a higher asymptote of the power spectrum at large wave 
numbers [Tj|, and in the halo approach, into a halo profile 
with more substructure and more mass in the center 
Note that the case described here is very similar to a 
comparison between an open and a flat model. 

This will show up in the weak lensing power spectrum 
in two ways. First, the earlier entrance in the non-linear 
regime, once projected, gives a transition between the 
linear and non linear regime at smaller t. Second, the 
amplitude of the power spectrum at small scale is ex- 
pected to be higher. BB01 proposed estimations of these 
two effects. In particular, due to the different evolution 
in the non-linear regime, the modification of the asymp- 
tote height is expected to go as the third power of the 
normalization modification in the linear regime. 

Probing the small scales of the shear power spectrum 
will reduce the degeneracy between the determination 
of dark energy properties and the normalization of the 
matter density fluctuation. The extend to which this can 
be done is the main subject of sec. ITTI 

Before going to next section, we would like to empha- 
size the fact that any modification in the nature of dark 
matter will be degenerated (at the level of the growth of 
structures of course) with the effect of dark energy. In 
particular, inclusion of hot dark matter will also modify 
the rate of growth. It is expected that this modifica- 
tion should decrease the amount of small scale structure, 
thus suppressing the effect of dark energy. We are then 
likely to underestimate the effect of dark energy in those 
models. 



E. Dark Energy Model 

The evolution of wq is a priori free. It has to be 
fixed by a proper model of dark energy. Several mod- 
els have been studied. The simplest being the minimal 
quintessence model, where wq is constant. Another very 



interesting class of models are the so-called tracking po- 
tential models. 

These models have been extensively described 
Their interest lie in the fact that their equation of state 
parameter is constant during most of the universe evo- 
lution. The constant equation of state is an attractor 
solution for wq when the expansion of the universe is 
dominated by another component (like the radiation or 
the matter). Of course, when the dark energy reaches 
the order of magnitude of the other energy densities, it 
leaves its attractor evolution. This attractor ensures that 
the initial conditions of dark energy does not have to be 
finely tuned. Whatever is the starting value of the dark 
energy 2 , it has to reach the attractor and will always 
exit the domination of matter at the same point, which 
in turn depends of the exact model. 

This explains the interest these models have met 
among the high energy physics community. In particular, 
it has been shown that some tracking potential models 
can be built within particle physics models. For exam- 
ple, P. Brax and J. Martin [2a| proposed a version of 
the Ratra-Peebles model [2U that can be embedded in 
super- gravity models. 

Minimal and tracking models are not the only dark 
energy models available. The problem however is that 
it does not exist any common framework to describe the 
different dark energy models that allows a direct compar- 
ison between them as we plan to do here, in the context 
of weak lensing. We attempt here to provide a suitable 
parameterization that will allow to constrain the dark 
energy properties using shear measurements. This is a 
simplification of the problem: we only have to model the 
evolution of dark energy that can leave an imprint on 
the weak lensing power spectrum. As a consequence, we 
only have to consider its impact on the growth of struc- 
ture and on the relation distance- redshift. 

We choose to parameterize the evolution of dark en- 
ergy in term of its equation of state. This choice is the 
most prevalent one. This is by no mean the only possi- 
ble parameterization |2~3] . As stated in section ITitl the 
knowledge of the EOS of dark energy is sufficient to solve 
Eq. <|112ll and to compute x( z ) an d D+- 

If we do not make any further assumptions, wq can 
freely vary between —1 and 1. The easiest solution is to 
assume that wq can be written as a power series of the 
redshift 

Wq = Wi z\ (23) 

This is the way most SNIa data are analyzed. More 
precisely, using the fact that the data collected relates 
directly to the integral of the expansion factor up to a 
redshift of order one, it is enough to only investigate this 



1 read fitted on N-body simulations 



2 Usually the initial condition are free within a few tens of order 
of magnitude. 
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equation of state in a perturbative development 

WQ = Wq + W\Z . . . 



(24) 



The possible determination of the two first orders of this 
development as Deen studied. 

This approach is not valid in our case. As noted above, 
one can only compute the growth of structure from re- 
combination. A perturbative development as Eq. 112411 is 
of course not suitable for our purposes; it leads to arbi- 
trarily growing equation of state, which is not physical! 
The full power series is also useless. As explained above, 
it is in the transition between linear and non-linear scales 
that one can expect to gain some information on the dark 
energy. As shown figure H and El one can expect some 
level of degeneracy even between a given varying and a 
non varying EOS models. Accordingly, it is very im- 
portant here that we greatly reduce the parameter space 
allowed to wq. 

Since we obviously cannot explore completely the equa- 
tion of state space if we want to restrict ourselves to a 
small number of parameters, we have to make some as- 
sumptions on the behavior of the dark energy equation 
of state. There are already too many different models of 
dark energy, and it is impossible to describe them with 
just one simple parameterization. Attempts have been 
made to generically describe dark energy with a simpler 
parameterization than the naive power series. They how- 
ever produce results with yet to many parameters j^| 
for our purpose. Another possible approach would be a 
principal component analysis designed for shear power 
spectrum inspired from ideas proposed in |23|. 

Here, to reduce the complexity, we add a physical as- 
sumption. We will only be interested in models who ex- 
hibit a behavior similar to the one of the tracking po- 
tential models. This is a very strong assumption. If the 
choice of model is hard to justify, we will however give 
here a few arguments in favor of the tracking potential 
behavior. 

The behavior of the dark energy equation of state at 
large redshift is in fact quite irrelevant for us. Indeed, 
tracking models assures that the equation of state of the 
dark energy is constant as long as it is dominated by the 
other components j^H- When dark energy is dominated 
its variation have little or no impact on the expansion 
of the universe. The growth of structure is not greatly 
modified by variation of wq when the dark energy com- 
ponent does not contribute significantly to the expansion. 
Accordingly, assuming a constant or varying equation of 
state during domination of matter or radiation is not rel- 
evant to our problem. Of course if one assumes that dark 
energy can be dominant at high redshift, this discussion is 
not valid. However, an such a dark energy model would 
leave a huge imprint on the CMB and would be most 
probably already ruled out by observations. 

At low redshift, when dark energy reaches the order of 
magnitude of the energy density of matter, it will start to 
contribute to the expansion, and induce a new period of 
acceleration. Variations of wq then leave a potentially 
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Figure 3: Comparison between an explicit SUGRA model and 
its parameterization. The equations of state are presented on 
the top panel, whereas the energy density, normalized to the 
critical density are on the bottom. Plain line is the SUGRA 
a — 6 model, long dashed is the logtan parameterization, 
short dashed, the log parameterization (see Eq. 1251 and 
l|26|l ). The parameter wq = —0.84, wi = 0.32 are measured 
on the SUGRA model. The log parameterization quickly fails 
to fit the equation of state above z ~ 1. It keeps a relatively 
good agreement on the dark energy density the to a higher 
redshift. It is not unexpected as the dominant contribution to 
H 2 is already the matter energy density. Thus slight variation 
on the equation of state of the dark energy are softened on 
SIq. The discrepancy, however build up very quickly to a 
factor 2 around z ~ 8. While not being in perfect agreement 
with the SUGRA model, the logtan parameterization does a 
better job at following Qq. 



strong imprint on the shear power spectrum, through 
modifications of the cosmological distances and structure 
growth. This is where our assumption on the shape of the 
model is important. Other models of the dark energy can 
have very different behavior. Our choice to only consider 
tracking like models theoretically reduce the reach of our 
conclusion, at least in its details. 

It has been shown j^JJ that the equation of state of 
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tracking models can be fitted at low redshift by a log 
function 



WQ 



Wo 



wi log(l + z). 



(25) 



This behavior is roughly valid up to redshift z ~ 1 at least 
for SUGRA and Ratra-Peebles potentials. Note that this 
equation of state parameter admits Eq. 112411 as its Taylor 
expansion at small z. 

The parameterization given in Eq. (J2SJ) fails quickly 
above z ~ 1 (see figure ) . It has to be modified to ac- 
count for the evolution of the equation of state for z > 1 . 
We propose to use an arctangent which has the property 
to quickly reach a constant value. We will then use 



1.003 - £> real 




(26) 



ivq + w\ log(l + z), if z < 1 
wo + wi [log (2) — arctan(l) - 



arctan(z)] , if z > 1 



This is a crude assumption, because the value of w at 
high z will be wo + wi(log(2) — arctan(l) + tt/2) which 
has no reason to correspond to the high z asymptote of a 
given tracking potential model. However, as said above, 
the dark energy component is completely dominated at 
high redshift by the energy density of matter, and the 
expansion is nearly insensitive to its evolution. In the 
end, the only thing that matters is the evolution of dark 
energy from the epoch when it starts to dominate. This 
epoch can be at redshift as high as 10. For example, 
in a SUGRA a — 6 model[lE| the energy density of dark 
energy represents 10% of the total energy density as early 
as redshift z ~ 5 (see figs. ElandQ]). 

The parameterization, Eq. j2Sj, is not very good at 
fitting the equation of state. As shown figure for the 
example of a SUGRA a = 6 model, it reasonably agrees 
with the dark energy density predicted in this model. 
This is not the relevant comparison anyway, as we should 
compare the growth and distances of the different models. 
This is done fig. 0] 

The agreement with the explicit tracking models we 
tested is around 3%. We will assume that this order 
of magnitude of the error holds for any other tracking 
model. 

We use this parameterization in our analysis. Although 
it is a ad hoc choice, it conveys most of the feature of a 
tracking model, and fits them, at least at the level of the 
growth of structure and the cosmological distances. It 
also has the advantage to be described by only two pa- 
rameters, wo and w\. We have fixed the change between 
the log and tan branches to z c = 1. Small variations 
of z c translate into small modifications of the growth of 
structure. For example, taking z c = 1.5 translates into 
one percent modification in D + /a, comparable with the 
error on the modeling. Of course, since we are only in- 
terested to lensing effect to redshift around 1, it leads 
to negligible modifications in the cosmological distances 



1 - \"+ 



Figure 4: Comparison between real tracking models and their 
approximated version using Eq. 12tifl . Top panel present the 
comoving distances, bottom the growth of structure. The 
lines plain line is a SUGRA a = 6 model, the dashed one, 
SUGRA a — 11 and the dotted a Ratra-Peebles, a — A model. 
The discrepancy on the angular distances computed with the 
real model and our parameterization is below the percent up 
to z = 4. The discrepancy for the linear growth if of order 3 % 
up to the recombination. Our approximated formula with its 
very small number of parameters gives a good approximation 
of the quantities on which is computed the weak lensing effect. 



and the projection effect 3 . Finally, since the parameter- 
ization 126t admits Eq. Il24t as its Taylor expansion, our 
results are directly comparable with the well advertised 
SNIa ones @. 



3 We do have a small dependency on redshift higher than 1 through 
the broad distribution of the source Ps(z). This effect is small 
enough to be neglected here. 
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II. RESULTS 

We perform a maximum likelihood analysis of the aper- 
ture mass statistic for a set of dark energy models. The 
method is well known and has been formerly described 
in Section fll Al describes the models and surveys 

that will be investigated. Finally, we present section III Bl 
the numerical results and a discussion on the degeneracy 
between the parameters. 



A. Parameter estimation 

We know from previous studies that the gravitational 
lensing by large scale structures depends mainly on four 
parameters: the matter energy density flo, the mass 
power spectrum normalisation as, its slope, and the red- 
shift of the sources As described above, we can 
safely ignore modification to the Cold Dark Matter trans- 
fert function due to dark energy. We thus use it and 
describe the slope of the power spectrum by the param- 
eter r. 

The lensing effect is also sensitive to other parame- 
ters, but to a lower extend. Therefore when studying 
the impact of the dark energy on cosmic shear, one has 
to incorporate the effect of the main parameters as well. 
We highlighted above that in the linear regime at least, 
dark energy variation should be degenerated with the 
normalization of the power spectrum. By keeping the 
main parameter we will be able to test other possible 
degeneracies. 

We assume that one of the main parameters is known. 
Indeed, the forthcoming lensing surveys are supposed to 
provide an accurate measurement of the distribution of 
the sources from photometric redshifts. Moreover, the 
redshift dependence is very similar to the as dependence. 
Hence, we will assume z s known. A slight error on its 
value can be translated in our result in broader as con- 
straint. 

Our set of free parameters is chosen as p = 
{wo,Wi,Q,Q,as,T). We deliberately choose a flat prior. 
The current CMB results are in very good agreement 
with a flat geometry 

We compute the likelihood £(p|d), where the data vec- 
tor d is the aperture mass (M|p) as function of scale: 



C = 



(2tj-)™/ 2 |S 



1/2 



exp 



(d-s) T S" 1 (d-s) 



(27) 



where s is the fiducial model vector and S := 
((d-s) (d-s)) is the covariance matrix. The covariance 
matrix is computed following the method described in 
[421 ] assuming the Gaussian field approximation. In this 
work we are interested in two surveys: one is the ground 
based Canada France Hawaii Telescope Legacy Survey 




0[arcmin] 

Figure 5: Aperture mass variance as function of scale for 
model3 (solid) and model2 (dashed), for the CFHTLS (thick 
bottom lines) and SNAP (thin top lines). The error bars 
show the statistical and sampling errors, assuming aGaussian 
statistic for the sampling error. 





deg 


n gal 


°e 


CFHTLS 0.9 


1790. 


20. 


0.44 


SNAP 1.2 


300. 


100. 


0.32 



Table I: Lensing surveys that will be part of the CFHTLS and 
SNAP projects (see text). Entries are source mean redshift 
z s , survey total area #j e g, source galaxy number density (per 

arc-min ), and intrinsic ellipticity dispersion ere- 



4 , and the other the spatial Super-Novae Acceleration 
Probe 5 . The observational properties of the lensing sur- 
vey associated with these two projects are summarized 
in table |U 

For the two surveys, we selected three fiducial models 
(with a cosmological constant fl\ = 0.7): 

• modell: p : = (-1,0,0.3,0.9,0.24) 

• model2: p 2 = (-0.8,0,0.3,0.9,0.24) 

• model3: p 3 = (-0.8,0.32,0.3,0.9,0.24) 

The first model is a pure cosmological constant case. Sec- 
ond is a minimal dark energy model with no variation of 
the equation of state. This type of model is widely used 
in the literature. From the discussion of section llDl it is 
expected that this kind of model under-evaluate greatly 
the effect of a varying EOS with identical final value. 
Last model has a varying EOS. The value of Wi has been 
choose so as to agree with a,a — 6 SUGRA model. It cor- 
respond to a strongly evolving equation of state model. 



4 http://www.cfht.Hawaii.edu/Science/CFHLS/ 
6 http://snap.lbl.gov/ 
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Models with a smaller w\ between our choices interpolate 
between our model2 and model3. 

We also have to make a choice on the range of param- 
eters we want to investigate. Maximum likelihood anal- 
ysis with five parameters is already a computationally 
expensive task. It can be reduced in part by narrowing 
the range of the parameters and the number of points in 
each direction. 

For the CFHTLS analysis, it is expected that we will 
mildly constrain the parameters. We thus used a rela- 
tively sparse grid and relatively wide parameter ranges. 
The matter density fly will be allowed to vary between 
0.1 and 0.5, while a s will be free between 0.6 and 1.1. 
The choices for this two parameters are quite conserva- 
tive. They allow to probe the full one sigma contour. The 
slope of the power spectrum is weakly constrained by the 
weak lensing measurement we probe its values between 
0.08 and 0.4. The results below ( figures I3ITT1) shows that 
this is more that enough to correctly probe the parameter 
space. 

For SNAP, we greatly reduce the range of parameters. 
The precision required here forces us to increase the num- 
ber of computed models, in particular in the f^o, &8 space. 
We thus suppose that it is enough to probe fio between 
0.28 and 0.32, and a 8 between 0.85 and 0.95. Never- 
theless, it is expected that by the time SNAP will collect 
data, previous weak lensing measurements, CMB, galaxy 
and cluster surveys will have cut down the accuracy on 
this parameters to these levels. We conservatively keep a 
relatively wide range on T (0.1 to 0.3, in agreement with 
the results for CFHTLS). 

For both models, we probed the dark energy param- 
eter space between -1 and 0.6 for wo and and 0.4 for 
w\. Note here that as described higher, we do not take 
into account models with an equation of state more neg- 
ative than -1. The upper bond on wq correspond roughly 
to the degeneracy expected between our target varying 
equation of state and a constant equation of state model 
(see section llD|l . We do not investigate negative wi mod- 
els. Negative Wi models are very close to the cosmological 
constant case, and should be strongly degenerated with 
it. It is very dubious that weak lensing will be able to dis- 
tinguish between them. The wi upper bond correspond 
to strongly varying equation of state. It is very difficult to 
reach this bond with SUGRA or Ratra-Peebles models. 



B. Numerical results — Discussion 

We first compute the aperture mass for our dark energy 
models. Figure presents the results for SNAP and the 
CFHTLS surveys. It shows that the evolution of the 
dark energy can lead to a 10 to 20% effect at small scale. 
As emphasized section ITTfl this is precisely the expected 
effect. 

Next we perform the likelihood analysis on our tar- 
get models. Figures IfllllOl and El show respectively the 
parameter predictions for the models 1, 2 and 3. All pos- 



sible combinations of pairs of parameters are plotted in 
order to show the direction of degeneracies. On each plot, 
the two hidden parameters are assumed to be perfectly 
known. We first note the strong degeneracy between the 
dark energy parameters (wq, wi) and the other param- 
eters. The full degeneracy between w\ and T is under- 
stood by the fact that the shape parameter describes the 
slope of the power spectrum, for a fixed normalisation 
08. Changing T will modify the ratio between linear and 
non-linear regime, and the scale of transition. As shown 
in section ITdI a change in w\ has a similar consequence. 

Even allowing for dark energy, the shear two points 
function remains a good constraint on 17 and a 8 . Figure 
shows the effect of unknown dark energy parameters 
(marginalised on wq and wi) on the measurement of fig 
and as- For the pure cosmological constant model (top 
panel), we see that the most probable models correspond 
to higher Qo and lower as than the fiducial model. For 
the model2 (bottom panel), the normalisation is under- 
estimated. This figure shows that the width of the fio, 
as contours is not dramatically affected, but the most 
probable models are changed. 

Super-Novae luminosity surveys have a small sensitiv- 
ity to the variation of the equation of state. In particu- 
lar, it is expected that without a strong prior on Qq they 
cannot provide much information on wi 0- 0. 01 ■ The 
question is whether the shear two points function also 
suffers from this kind of limitation or not. Figure Hand 
IHlshow the predictions for wq and wt, respectively for the 
CFHLS and SNAP observations. The left panels corre- 
spond to the contours obtained with a perfect knowledge 
of Qo, as and T. The middle panels are for a known T, 
but marginalised over f2o and a%. The right panels are 
for known Slo and as, and marginalised over T. The top 
panels are for model 2, and the bottom panel for model 
3. The important result here is that the marginalisation 
over Qm and as do not increase too much the width of 
the contours, it only restores a degeneracy between Wq 
and w\. 

Contrary to the angular diameter distance tests, the 
weak lensing is sensitive to the evolution parameter w\. 
The marginalisation over T restores the degeneracy along 
a different direction, but still does not increase the con- 
tours width significantly. It means that even with a lim- 
ited knowledge on external important parameters, it is 
still possible to constrain the quintessence, in particu- 
lar when it evolves with time. In that case indeed (i.e. 
wi 7^ 0) , the increase of the lensing signal is large enough 
to allow the CFHTLS observations to rule out a pure cos- 
mological constant case. However, an accurate joint mea- 
surement of the quintessence parameters and the others 
is not possible using the lensing power spectrum alone, 
because of the strong degeneracy between Q and as ■ 

This degeneracy is broken with the SNAP lensing sur- 
vey: according to Figure |H1 one sees that cosmic shear 
observations alone with the SNAP satellite, provide con- 
straints which are competitive with the SNIa constraints 
from the same satellite. The expected constraints from 
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Figure 6: Contours in the do, ""8 space when marginalised 
over the quintessence too € [—1, —0.7] and toi € [0,0.4]. This 
is given for the CFHTLS experiment, left panel corresponds 
to model 1, and right panel to model 2. 



SNIa alone, assuming a perfect knowledge of £lo, is 
sketched on this figure (solid line) . It shows that SNIa are 
less sensitive to w\ than weak lensing. Therefore a com- 
bination of SNIa and cosmic shear could simultaneously 
probe the dark energy and its evolution. More precisely, 
figures EH 1101 and El show that the knowledge of wq is 
irrelevant for constraining Slo from cosmic shear. On the 
other hand, the SNIa measurements are degenerate be- 
tween wq and flo ■ A combination of the two experiments 
provide a simultaneous measure of wq and Slo without 
the need for an external measurement of Slo- We can 
then use the lensing constraints on wo and wi (Figures 
[3 and [HJl to estimate the dark energy evolution w\. In 
fact even a poor knowledge of f^o can be tolerated; we 
known from Q that a marginalisation over of the SNIa 
measurements increases the wq, wi contours perpendic- 
ularly to the increase of the same contours from cosmic 
shear with poor knowledge on flo (Figures and |H1 mid- 
dle panels). Adding the cosmic microwave background 
over-constrains the parameter space, because the con- 
tours in the flo, wq space are 'perpendicular' to the SNIa 



Figure 7: CFHTLS constraints with lensing alone on too and 
toi. Top panels: model2, bottom panels: model3. Left plot is 
assuming all other parameters are known (see Figure flOl and 
llljl . Middle plots is when the mean density and the power 
normalisation are marginalised (flat prior) over fio £ [0.1, 0.5] 
and as £ [0.6, 1.1]. The right plots show the contour for the 
marginalisation F £ [0.1,0.4]. 




Figure 8: Same as Figure [Jfor the SNAP survey. The 
marginalisation is performed over the intervals Slo £ 
[0.28,0.32] and cr 8 £ [0.85,0.95] for the middle plots and 
r £ [0.1,0.3] for the right plots. The line show the direction 
of degeneracy of the SNIa (with the supposition of a perfect 
ilo knowledge) 



and cosmic shear constraints [4J|. Weak lensing, cos- 
mic microwave background and SNIa provide therefore 
an ideal set of complementary experiments for constrain- 
ing the dark energy beyond the constant energy density 
case because weak lensing measurement breaks the 
degeneracy with the dark energy evolution. 

Earlier work has shown that cosmic shear provides also 
independent constraints on Qq from the measurement of 
high order statistics of the convergence 
Dark energy modifies mildly this picture. At the level of 
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the quasi-linear regime, it only affects the three points 
function of the convergence field through the projection 
effect 0|- The modifications are expected to be more 
important at small scales j^J]. One can see from Figure 
131 1101 and El that this additional information is not nec- 
essary, given the degeneracy among wq in the (wq, Sl ) 
space. However, such external constraint could be very 
helpful to pin down the degeneracy with erg, and con- 
sequently to reduce the degeneracy between Wq and w\, 
helping to narrow the constraint on Wi . 

III. CONCLUSION 

We investigated the possibility to constrain the evolu- 
tion of dark energy evolution from measurements of the 
gravitational lensing by large scale structures. We used 
the fact that the non-linear growth rate of structures is 
significantly affected. This is a consequence of an earlier 
influence of dark energy on the expansion of the Uni- 
verse. It was found that the cosmic shear effect is an 
ideal probe of the evolution of dark energy, in opposition 
with experiments based on angular diameter distances 
like SNIa and cosmic microwave background, which are 
better suited to measure the "constant" part of dark en- 
ergy equation of state (in a particular parametrization) . 
The degeneracy with other parameters (O , o"8 and T) 
restore a degeneracy between u>q and w\, but the width 
of the contours in that space are slightly affected. There- 
fore a linear combination of wo and wi is well measured 
using weak lensing, even with a poor knowledge on Qq 
and ag. 

It is generally believed that the measurement of the 
dark energy equation of state parameter as a constant is 
such a difficult task, that we should not even dream to 
measure its evolution. We have shown here, for a class 
of models, that the sensitivity of the large scale struc- 
tures to a simple evolution parameter w\ is as easy (or 



difficult!) as wq to measure. Consequently, we found 
out that a combination of cosmic shear, SNIa and cosmic 
microwave background provide orthogonal constraints of 
the parameters wq, w\ and f^o, which opens great oppor- 
tunities to probe non-trivial models of dark energy. For 
the set of models studied here, we found that these three 
experiments over constrain these parameters. 

One should note that the difference in the amplitude 
of the cosmic shear signal between model 2 (wq = —0.8, 
wi = 0) and model 3 (w = —0.8, Wi = 0.32), at scales 
below 5' reaches 10%. This is large compared to the sta- 
tistical errors of the CFHTLS and SNAP surveys. How- 
ever, it is yet within the limits of the Point Spread Func- 
tion (PSF) correction and non-linear modeling accuracies 
[111 ]. If one wants to measure the dark energy evolution 
as proposed here, it is clear that we need to perform ray- 
tracing simulations for the class of models we want to 
investigate, in order to calibrate the non-linear modeling 
[filj . The PSF correction is an entirely different issue, 
which is not discussed here, but there is good hopes to 
be able to reduce the systematics level by a factor of 5 
to 10 t 52|, which should be enough for our purpose here. 

Redshift degeneracy was not discussed, but it is not 
different from the a$ and the flo degeneracies. What has 
been said for these parameters also applies to the source 
redshift. In the future, photometric redshifts will provide 
accurate source redshift measurements, as we do not need 
an accurate redshift for each lensed galaxy. 
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Figure 9: Constraints obtained with the CFHTLS survey for a cosmological constant model (modell). We assumed strong 
prior for the hidden (not shown) parameter on each plot. The cross represents the fiducial model. 
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Figure 10: Same as Figure Ejfor a quintessence, non evolving, model2 (see table 
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Figure 11: Same as Figure E3 for an evolving dark energy model (model3). 
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